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-L ■ Abstract. The blackbody spectrum of CMB was created in the blackbody photosphere at redshifts 

z > 2 x 10 6 . At these early times, the Universe was dense and hot enough that complete thermal 
equilibrium between baryonic matter (electrons and ions) and photons could be established on time 
scales much shorter than the age of the Universe. Any perturbation away from the blackbody spec- 
trum was suppressed exponentially. New physics, for example annihilation and decay of dark matter, 
can add energy and photons to CMB at redshifts z > 10 5 and result in a Bose-Einstein spectrum with 
a non-zero chemical potential (p). Precise evolution of the CMB spectrum around the critical redshift 
of z - 2x 10 6 is required in order to calculate the /i-type spectral distortion and constrain the underly- 
£X ing new physics. Although numerical calculation of important processes involved (double Compton 

-r^ ■ process, comptonization and bremsstrahlung) is not difficult with present day computers, analytic 

solutions are much faster and easier to calculate and provide valuable physical insights. We provide 
precise (better than 1%) analytic solutions for the decay of /i, created at an earlier epoch, including 
all three processes, double Compton, Compton scattering on thermal electrons and bremsstrahlung in 
the limit of small distortions. This is a significant improvement over the existing solutions with ac- 
/^ ', curacy ~ 10% or worse. We also give a census of important sources of energy injection into CMB in 

standard cosmology. In particular, calculations of distortions from electron-positron annihilation and 
primordial nucleosynthesis illustrate in a dramatic way the strength of the equilibrium restoring pro- 
cesses in the early Universe. Finally, we point out the triple degeneracy in standard cosmology, i.e., 
the fi and y distortions from adiabatic cooling of baryons and electrons, Silk damping and annihilation 
of thermally produced WIMP dark matter are of similar order of magnitude (~ 10~ 8 - 10 -10 ). 
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1 Introduction 

There are several important events in the history of the Universe which provide the foundations for 
the standard cosmological model. The first event is the big bang, which created the present expanding 
Universe filled with matter, radiation and dark energy. One of the goals of cosmology is to reconstruct 
these initial conditions and learn about the high energy physics at early times. We should emphasize 
here that the blackbody spectrum of cosmic microwave background (CMB) is not an initial condition. 
The blackbody spectrum is created and maintained dynamically throughout the early histoiy of the 
Universe (z > 2 x 10 6 ) by standard physics processes and how this happens is the main topic of the 
present paper. The CMB spectrum thus provides information about the physics of the Universe after 
the big bang. 

The other important events in standard cosmology are: the formation of light elements in pri- 
mordial nucleosynthesis [1], recombination of electrons and ions to neutral atoms [2, 3], and reion- 
ization of the Universe [4] by the radiation emitted by first stars and galaxies. We should also mention 
electron-positron annihilation at 10 10 > z > 10 8 which more than doubles the entropy of CMB and 
raises its temperature by ~ 40%. The important events in the history of the Universe are sketched 
in Fig. 1. Big bang nucleosynthesis (BBN) theory together with measurement of light elements 
abundances constrains the photon to baryon number density at z ~ 10 8 [see 5, for a review] and is 
the earliest direct evidence and measurement of electromagnetic radiation in the early Universe. The 
fact that the photon energy density inferred from BBN and Cosmic microwave background (CMB) 
at z ~ 10 3 [6] are close to each other (consistent within ~ 2 - cr) also implies that we do not have 
arbitrary freedom in adding energy to CMB between these two epochs. There is, however, no di- 
rect way to constrain the energy density in photons before BBN and the epoch of electron-positron 
annihilation. 
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Figure 1. Important events in the history of the CMB spectrum and anisotropy formation in big bang cos- 
mology. Redshift range (2 x 10 6 > z > 10 5 ), where the energy injection would give rise to a Bose-Einstein 
spectrum (yU-type distortion), is marked as fi. At much smaller redshifts (z < 10 4 ), any heating of CMB through 
Compton scattering would create a y-type distortion. The spectrum in the intermediate redshift range would 
not be a pure /u or y type but in between the two types. 



The CMB (anisotropy and polarization) is at present the most precise cosmological probe. The 
CMB spectrum was created at z > 2 x 10 6 , and this critical redshift defines the blackbody surface for 
our Universe. Spatial fluctuations in the temperature of the CMB were imprinted much later, when the 
electron and protons recombined to form hydrogen atoms at z ~ 1100. This second boundary defines 
the well known last scattering surface, the structure of which is encoded in the photon visibility 
function [first studied by 7]. The anisotropics and physics at the last scattering surface have been 
very well studied and accurate analytic [8] and numerical solutions [9, 10] have been available for 
some time, motivated by the precise experiments such as WMAP [6] and Planck [11]. 

The blackbody spectrum, once created at high redshifts (for example before the time of electron- 
positron annihilation), is preserved by the adiabatic expansion of the Universe at all subsequent times. 
However, if there is energy release at lower redshifts , for example by particle decay and annihilation 
or Silk damping, it will distort the CMB spectrum away from the Planck form. In this case, Zeldovich 
and Sunyaev [12] first demonstrated that bremsstrahlung alone cannot recreate blackbody spectrum 
until very high redshifts, almost up to the time of electron-positron annihilation. The problem of evo- 
lution of the CMB spectrum through the blackbody surface, in the presence of heating, was solved 
analytically by [13] including the processes of comptonization and emission and absorption of pho- 
tons with special emphasis on bremsstrahlung. Comptonization is the process of redistribution of 
photons over frequency, resulting from the Doppler and recoil effects of Compton scattering of pho- 
tons on thermal electrons [14]. Since the double Compton cross section [15, 16] has a dependence 
on frequency similar to bremsstrahlung, the solution of Sunyaev and Zeldovich [13] also allowed 
inclusion of double Compton emission and absorption, which is dominant over bremsstrahlung in a 
low baryon density Universe such as ours and was first considered by [17]. Double Compton emis- 
sion or absoiption is just the first radiative correction to the process of Compton scattering just as 



bremsstrahlung is the first radiative correction to the scattering of electrons on nuclei [see discussion 
in 18]. The approximations used in these analytic solutions result in accuracy of ~ 5 - 10% at 
redshifts (z < 3 x 10 6 ) and worse than 10% at higher redshifts. Much better accuracy is, of course, 
achievable in numerical solutions [19-23]. 

The type of spectrum we get for the CMB is determined by the Compton y parameter, defined 
in Eq. (2.5). Bose-Einstein spectrum can be created only at high redshifts z > I0 5 ,y > 1 [7], 
when comptonization of CMB is very efficient. Any addition of energy and photons to CMB at 
z < 10 5 ,y << 1 inevitably distorts the CMB spectrum with a v-type distortion [12]. The v-type 
distortions, in contrast to yu-type distortions, can be created at low redshifts up to z - 0. Reionization 
at z ~ 10 in particular is expected to create y-type distortions of magnitude y e ~ T r ikQT e /(m e c 2 ) ~ 
10~ 7 where r n - « 0.1 is the optical depth due to reionization, T e « 10 4 is the average temperature 
during reionization, k% is the Boltzmann's constant, c is the speed of light and m e is the mass of 
electron. Due to uncertainties in reionization physics, it will be extremely difficult to separate these 
low redshift distortions from the .y-type distortions created before recombination. COBE FIRAS 
experiment [24] measured the CMB spectrum with high precision and placed a constraint on the 
chemical potential of Bose-Einstein spectrum for CMB of /i < 9 x 10 5 , thus confirming that the 
CMB has a Planck spectrum at very high accuracy. 

In standard model of cosmology we can get constraints on the energy density of radiation from 
two distinct and very precise observables. The first is the deuterium abundance, which gives the 
baryon number to photon number ratio rj = (5.7 ± 0.3) x 10 -10 during primordial nucleosynthesis at 
4 x 10 8 > z > 4x 10 7 [5, 25]. The second is the measurement of CMB anisotropics, which constrains 
the baryon to photon ratio 77 = (6.18 ± 0.15) X 10~ 10 [6] during recombination at z ~ 1 100. The fact 
that these two values of baryon to photon ratio are almost identical with small error bars means that 
we do not have arbitrary freedom in adding energy to CMB, for example, with the introduction of 
new physics. In fact any addition of energy/entropy to CMB between primordial nucleosynthesis and 
recombination cannot be more than a small percentage (< 10%) of the photon energy density during 
BBN. 

In standard cosmology, the chemical potential of CMB is expected to have a magnitude of 
// ~ 10~ 9 - 10~ 8 resulting from the heating of CMB by Silk damping. A similar magnitude but 
opposite sign distortion, /1 « -2.7 x 10 9 , is expected from cooling of photons due to energy losses to 
baryons and electrons, which have a different adiabatic index (5/3) than radiation (4/3) and as a result 
cool faster than the radiation as the Universe expands [23, 26]. A detection of a chemical potential of 
magnitude greater than 10~ 8 will therefore mean existence of non-standard physics (or a small scale 
power spectrum that is bigger than what is expected from extrapolating the large scale power spectrum 
as measured by CMB and large scale structure) at these redshifts. Proposed experiment Pixie [27] has 
exactly this level of sensitivity. Constraining high energy physics using yu-type distortions requires 
precise calculation of evolution of the CMB spectrum through the blackbody surface at z ~ 2 X 10 6 . 
Analytic solutions provide valuable physical insight in addition to being much easier to compute. 
Motivated by these factors we try to improve the method of [13] to achieve better than 1% precision 
in analytic solutions. These solutions should prove useful in predicting the signal from models of 
high energy physics which can provide energy injection to CMB at these high redshifts. 

In the last part of the paper we apply our analytic solution to examples from the standard model 
of cosmology including electron-positron annihilation and dark matter annihilation. These simple 
examples in particular demonstrate the efficiency of processes responsible for maintaining thermal 
equilibrium in the early Universe. They also demonstrate the difficulty of creating spectral distortion 
in CMB at high redshifts. There remains thus only a narrow window at redshifts 10 5 < z < 2 x 10 6 
when a Bose-Einstein spectrum can be created. 



We use WMAP [6] best fit parameters for ACDM cosmology for numerical calculations. 



2 Thermalization of CMB 



Compton scattering is responsible for creating a Bose-Einstein spectrum of photons if the rate of 
comptonization (i.e. the redistribution of photons over the entire spectrum by Compton scattering) is 
greater than the expansion (Hubble) rate of the Universe [12]. This condition is satisfied at redshifts 
z > 10 5 . In the non-relativistic regime comptonization is described by Kompaneets equation [14]. 
Compton scattering, however, conserves photon number and therefore the spectrum obtained as a 
result of comptonization will in general have a non-zero chemical potential. At redshifts z > 2 x 10 6 
double Compton scattering, and to a lesser extent bremsstrahlung, can emit/absorb photons at low 
frequencies very efficiently because the optical depth and the absorption coefficient of these two 
processes increase with decreasing frequency as <x v~ 2 . Compton scattering is then able to redistribute 
these photons over the entire spectrum. The net result is that the chemical potential is exponentially 
suppressed and a Planck spectrum is established. 

The equilibrium electron temperature in a radiation field with occupation number n(x), where 
x = hvjkTy and T y is the temperature of reference blackbody, is given by [28, 29] 



J e J(n + n 2 )x 4 dx 
T y 4jnx 3 dx 



(2.1) 



The rate (ScrjEy/(3m e c) x n e /(n e + «,•)) at which electron/baryon plasma achieves equilibrium tem- 
perature T e given by Eq. (2.1) is shown in Fig. 2 in the topmost curve. Rates of bremsstrahlung 
absorption K\, r (e x,: - l)/x e 3 , double Compton absoiption K d c(e Xe - l)/*e 3 and Compton scattering 
Kq are also compared with the Hubble rate in Fig. 2. n e is the electron number density, aj is the 
Thomson cross section, T e is the electron temperature, «,• is the number density of ions, x e = hv/k^T^ 
is the dimensionless frequency corresponding to frequency v, h is the Planck's constant, E y is the 
energy density of photons. T e = T y = 7cmb(1 + z), 7cmb - 2.725K and blackbody spectrum is 
assumed for this figure. The rate coefficients Kq, ^dc> ^br are given by: 
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af s is the fine structure constant and «b is the baryon number density. I^c - J dx e x s 4 n(l + n) ~ 
25.976 for a blackbody spectrum at temperature T e , Q. b is the baryon density parameter, ho is the 



Hubble parameter and Yu e is the primordial helium mass fraction, g\, r - TiiZ£nigB(Zi,x e ,T e )/nB 
is the average gaunt factor for bremsstrahlung, «,• is the number density of ion species i and Z, the 
charge of ion. The sum is over all ionic species, which for the primordial plasma at high redshifts 
consists of protons and helium nuclei. Accurate fitting formulas for g^ have been provided by [30]. 
gdC is the gaunt factor for double Compton scattering and accurate fitting formula for it has been 
calculated recently by [23]. We use these fitting formulae for numerical solution. For reference at 
x e = 0.01, z = 2 x 10 6 , gdc = 1-005 and gbr = 2.99 and we use these values for our analytic solutions. 
These gaunt factors are slowly varying functions of time and frequency and can be assumed to be 
constant in the redshift range of interest for analytic calculations. We justify this assumption below. 
The division into j-type and Bose-Einstein regions depends on the Compton y parameter, 

y{z) = dz — ■ (2.5) 

Jo m e c H(l+z) 

During radiation domination (z » z e q = 3.2 x 10 4 ) the integral can be carried out analytically, giving 
y(z) » 4.9x 10~ n (l +z) 2 , where z e q is the redshift when matter energy density equals radiation energy 
density and we have assumed 3.046 effective species of massless neutrinos [31]. We have, for WMAP 
ACDM cosmological parameters [6], y(1.5 x 10 4 ) « 0.01, y{A.l x 10 4 ) « 0.1 andy(1.5 x 10 5 ) « 1. 
For y < 0.01 « 1 we have ay-type distortion and for y > 1 a Bose-Einstein spectrum. The spectrum 
attained for 0.01 < y < 1 is in between a pure .y-type and Bose-Einstein. We choose z = 5x 10 4 as an 
approximate division between the two types of distortions for our estimates. This division is accurate 
if the energy injection between 1.5 x 10 4 < z < 1.5 x 10 5 is small compared to the total energy 
injection at z < 2 x 10 6 . If most of the energy injection happens between 1.5xl0 4 < z < 1.5xl0 5 , 
then numerical calculations must be performed to get accurate final spectrum. We note that for small 
distortions, we can calculate distortions arising from different physical processes separately and add 
them linearly. 

2.1 Kinetic equation 

We will follow [13] in solving the kinetic equation for photon distribution using stationarity assump- 
tion and then do an iteration to relax this assumption to arrive at a more accurate solution. The kinetic 
equation for the evolution of photon occupation number n(x e ) in the presence of Compton scattering, 
double Compton scattering and bremsstrahlung is given by [see also 20, 21] 
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The first term with coefficient #c is the Kompaneets term describing Compton scattering. The three 
terms in the square brackets in Kompaneets term describe photon diffusion in frequency due to the 
Doppler effect, electron recoil and induced recoil effects respectively. The second term represents 
emission and absorption of photons due to double Compton (ic'dc) and bremsstrahlung (^br)- The 
last term in Eq. (2.6) arises because we are evaluating the time derivative at constant x e = hv/k^T c 
instead of constant frequency v, and electron temperature changes with time 1 because the photon 
distribution is evolving (see Eq. (2.1)). 

We are interested in the regime where comptonization is efficient and deviation from a Planck 
spectrum is small. In this regime the photon spectrum is described by a Bose-Einstein distribution 



'This is in addition to the usual (1 + z) dependence due to the expansion of the Universe. The variable x c is invariant 
w.r.t. the expansion of the Universe. 
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Figure 2. Rates of bremsstrahlung, comptonization and double Compton scattering. Bremsstrahlung and 
double Compton absorption have a dependence on frequency ~ x~ 2 . Rates above are for x e - 0.01 which is 
the typical value of critical frequency at which double Compton/ bremsstrahlung absorption becomes equal to 
comptonization rate (Eq. (3.2) and Fig. 3). Also shown is the rate at which electron/baryon plasma achieves 
equilibrium temperature T e given by Eq. (2.1). T e - T y and blackbody spectrum is assumed for this figure. 
Compton scattering is able to create Bose-Einstein spectrum at z > 10 5 while for redshifts z < 15000 any 
energy injection gives rise to a v-type distortion [12]. The distortion has a shape in-between y-type and //-type 
for redshifts 1 .5 x 10 4 < z < 10 5 and we choose z - 5 x 10 4 to divide our estimates between /U-type and y-type. 



with a chemical potential much smaller than unity in magnitude. 

1 1 « 



n(x e ) = 



e x c +)i _ 1 g*e _ 1 ^( e v e _ 1)2 

The total energy density and number density of photons is then given by 2 
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where cir - n 3 p is the radiation constant, b R = - 3 B 3 ^ , f is the Riemann zeta function with 
f(3) ~ 1.20206, h = Jx 3 n pl (x)dx = n 4 /\5,I 2 - J x 2 n pl (x)dx = 2^(3), and n pl (x) = l/(e x - 1). 



2 The chemical potential is not constant but a function of frequency at low frequencies but this dependence is important 
only for calculating the photon production rate. We can ignore the frequency dependence in calculating the total energy 
and number density in the spectrum since the contribution from low frequencies (x <K 1) to these quantities is small and 
the constant ji assumption introduces negligible error. 
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In order to cancel the effect of the expansion of the Universe we will use the blackbody spectrum 
with temperature T y = 7cmb(1 + z), E y is its energy density and N y its number density. If we have 
a source injecting energy density at a rate £ and photon number density at rate N, where £ = E/E y , 
N = N/N y , we get using Eqs. (2.8) and (2.9) and with |//| « 1 

±JL) s s =4 ±Jh)-sm«ji (2 , 0) 

dt \Ey) £ At \T y ) h dt 

dt \N y ) N dt \Ty) 3/ 2 At v ' 

For small distortions E « E y ,N « N r and we have at lowest order, £/£ « £ and Af/Af ~ N. The 
photon production due to double Compton and bremsstrahlung can be calculated by taking the time 
derivative of Eq. (2.9). Multiplying the kinetic equation Eq. (2.6) by x e 2 and integrating over x e and 
using it in the the time derivative of Eq. (2.9), the terms involving dT e /dt cancel out and only the 
bremsstrahlung and double Compton terms contribute, giving 

d (N\ 1 f°° e" Xe 

tM^T ) = T dx e (K dC + K hr ) 1- B (e*-1)] 

At \N y ) I 2 Jo x e 

(K dC + K br ) r n 

« -j d *e— — -■ (2.12) 

h J x e (a* - 1) 

In general there may be additional sources of photons, for example, if energy is injected as an elec- 
tromagnetic shower resulting from decay of a heavy particle, the resulting cascade may produce 
non-negligible amount of photons. However in most cases of interest and at high redshifts this is 
much smaller than the photon production from bremsstrahlung and double Compton. 

2.2 Numerical solution 

In order to calculate the precision of the analytic formulae derived below, we compare them with 
the numerical solution of the coupled system of Eqs. (2.6) and (2.1). The initial spectrum for the 
numerical solution is a fi type distortion of magnitude 10~ 5 at z - 5 x 10 6 with the chemical potential 
decaying exponentially with decreasing frequencies at x « 1 (see Eq. (3.2) below). The results 
are not sensitive to the exact form of initial spectrum if the starting redshift is » 10 6 , what matters 
is the total energy input into an initial blackbody. We solve equation (2.6) iteratively in small time 
steps (using Compton y parameter as the time variable) of 6y = 0.1. In the first iteration we use the 
analytic solution for the evolution of electron temperature. In the second iteration we use the solution 
of first iteration to calculate the electron temperature using Eq. (2.1). As a frequency variable we 
use a x = hv/kftT where T is the reference temperature which evolves just by redshifting due to 
cosmological expansion and is equal to the electron temperature at the start of each iteration step. 
This gets rid of the last term in Eq. (2.6) and introduces factors of T e /T in the Kompaneets and 
bremsstrahlung/double Compton terms. We also write the total spectrum as a sum of blackbody part 
at reference temperature T and a distortion part, keep only the terms linear in distortions and solve the 
linearized Eq. (2.6) for distortions, as the zeroth order blackbody part vanishes. Implicit backward 
differentiation method is used to solve the PDE. We use logarithm of x as the second independent 
variable with the variable step size in the x direction. 

The main source of error is the deviation of electron temperature used in solving the PDE from 
the correct temperature given by Eq. (2.1), resulting in violation of energy conservation. In our 
solution, the maximum error in energy conservation as a fraction of energy in the distortion occurs 
at high redshifts, when there is strong evolution of ju, but even this is ~ 10 -5 in each iterative step. 



The error with respect to the total energy density in photons is, therefore, ~ 10~ 10 , since we start 
with an initial distortion of ~ 10 -5 , and is much smaller in the later steps as the distortion decays 
exponentially. An important point to note here is that since we change the reference temperature to 
the current electron temperature at the beginning of each iterative step, and solve and track only the 
distortion part, the error in // in individual iterative steps does not accumulate but is in fact suppressed 
in the subsequent evolution of the spectrum by the visibility factor (defined below). The solution 
obtained can thus be considered almost exact for the purpose of the present paper. 

The high frequency spectrum is forced to be Wien at 120 > x > 100 with the chemical potential 
and temperature given by analytic solution. The low frequency boundary is at x = 10~ 5 and the 
spectrum at x < 2 x 10~ 5 is forced to be blackbody with the temperature equal to the analytic elec- 
tron temperature. Since our boundaries are far away in the distant Wien/Rayleigh- Jeans tails, where 
there are negligible amount of photons/energy , the solution is not sensitive to the exact boundary 
conditions. Since our boundary conditions are approximately equal to the true solution, we are able 
to use large steps in the x direction resulting in considerable speedup in the numerical calculation 
compared to a calculation with arbitrary (but reasonable) boundary conditions. The initial spectrum 
is evolved until recombination (z = 1 100) although the // distortion at x > 0.1 is effectively frozen-in 
after z ~ 10 4 . Further details on numerical issues can be found in [20, 22, 23]. 

3 Solution using stationarity approximation 

In order to calculate the integral in Eq. (2.12) we need to know the occupation number n(x e ). Note 
that we cannot use Bose-Einstein distribution as an approximation since the integral in this case 
diverges at small x e . The reason is that double Compton and bremsstrahlung rates diverge at small 
frequencies and establish Planck spectrum. We can get an approximate solution from the kinetic 
equation Eq. (2.6) [14] by assuming that the instantaneous spectrum is stationary. Thus, neglecting 
time derivatives, making chemical potential function of frequency, 3 fi(x e ), and also assuming, x e « 
l,fi(x e ) « 1 Eq. (2.6) simplifies considerably, 

1 d 9 du u 

= -K c ——x & 1 -f- + {K dC + K hr )-^- A (3.1) 

The solution of this ordinary differential equation with the boundary condition /i(0) = is given by 
[13] 
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where fi c is normalization specified by chemical potential at large x e . Thus fi(x e ) decays exponen- 
tially at small frequencies and goes to constant at large frequencies. x c is also the frequency at which 
comptonization rate is equal to photon absoiption rate due to double Compton and bremsstrahlung. 



3 We note that any spectrum can be described by a frequency dependent chemical potential. 



Similarly we can also define a frequency at which the photon absoiption rate is equal to the Hubble 
rate H(z) [21], 
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;c c and x^ are plotted in Fig. 3. We note that x c « 1 in the redshift range of interest. This is 
consistent with assumptions made in the derivation. Also, since x c ~ 0.01 in the redshift range of 
interest and gaunt factors are slowly varying functions of frequency and temperature, we can assume 
gdc * gdc(O.Ol) = 1.005 and g br * 2.99. 

We should emphasize that the sole purpose of finding an accurate solution for the spectrum at 
low frequencies is to calculate precisely the photon emission/absorption due to bremsstrahlung and 
double Compton. In particular, only the spectrum at x e > 0.1 is frozen-in at x «; 10 5 . The low 
frequency spectrum continues to be affected by bremsstrahlung (at smaller redshifts, during recombi- 
nation and after), which tries to bring the spectrum in equilibrium with the electrons, which are cooler 
than the radiation due to adiabatic expansion [23, 26]. Furthermore, most of the photons are created 
around the critical frequency x c . The double Compton gaunt factor at the critical frequency gdc(-Xc) 
deviates from the value 1.005 by less than 0.5% in the interesting redshift range of 10 5 < z < 10 7 . 
The bremsstrahlung gaunt factor has a maximum deviation of ~ 7% from 2.99 in the same redshift 
range, but since it is only a small correction to the dominant double Compton process, the error in the 
final solution for /i evolution is small. Thus the assumption of constant gaunt factor is an excellent 
one for the present problem and is further justified by a comparison of the analytic and numerical 
solutions. 

We can now use the solution Eq. (3.2) to evaluate the integral Eq. (2.12) (ignoring the x e 
dependence of gaunt factors), 
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Equations (2.10) and (2.11) along with the above solution give the following equation for the evolu- 
tion of chemical potential p - p c at x > 1, 
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where C = 0.7768, B - 1.803. jr are the extra photons injected from processes other than the 

Extra 

low frequency bremsstrahlung and double Compton photons calculated above, for example, from the 
same source which injects energy. We assume this extra term to be negligible in the rest of the paper. 
The solution of Eq. (3.5) at z - is given by 
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where z m i n « 5 x 10 4 , z.i is the initial/maximum energy injection redshift and we have defined the 
effective blackbody optical depth (which is frequency independent at x > x c ) 
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Figure 3. The frequency x c and xh at which the Compton and Hubble rates respectively are equal to the 
sum of the bremsstrahlung and double Compton rate. If Compton rate is also greater than the Hubble rate, then 
bremsstrahlung and double Compton can establish complete thermodynamic equilibrium (blackbody spectrum) 
below x c otherwise complete thermodynamic equilibrium is established below xu- Above x c , at redshifts 
Z > z c ~ 10 5 , a frequency dependent chemical potential /i is established. At x » x c , the chemical potential 
has an almost constant (frequency independent) value jj c and we have thus a Bose-Einstein spectrum. The 
chemical potential decreases with time due to the photons created by bremsstrahlung and double Compton at 
low frequencies and redistributed by Compton scattering over the entire spectrum. 
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O r is the total radiation density parameter, and H(0) is the Hubble constant today. It is interesting to 
note that in the absence of double Compton scattering we would have z.br ~ 6 x 10 6 . The presence 
of double Compton increases the critical frequency x c from its bremsstrahlung only value, and thus 
reducing the bremsstrahlung emission. It is also straightforward to calculate, if needed, the chemical 
potential at any intermediate redshift z' > Zmin using Eq. (3.6) by replacing T{z) (and similarly for 
T{zi)) with T{z) - T(z') and also replacing the lower limit z m in in the integral with z' , 

Hiz) = Kzde-^- T ^ + CB f -JJL- (| - 1 £ ) e -[TO-™l (3.9) 

J z/ (l+z)H\B 3 N Ex J 

The solution given in Eq. (3.2) corresponds to Eq. (15) in [13] but including the double Comp- 
ton process. Similarly, the solution in Eq. (3.7) generalizes Eq. (20) of [13]. The dominant term 
(Zdc)in Eq. (3.7) is due to the double Compton process with the bremsstrahlung term (zbr) providing 
a small but important correction. 

4 Improved solution by approximating non-stationarity using previous solution 

We will see below that the solution arrived at in the previous section underestimates the photon 
production. It turns out that the stationary solution, which is normalized at high frequencies, un- 
derestimates the chemical potential at small frequencies where most of the photons are being pro- 
duced/absorbed. We find below the correction for the normalization of chemical potential, which 
enables us to improve the formula for blackbody optical depth Eq. (3.7) in the solution Eq. (3.6). 
The result of the computations, using the new analytic formula Eq. (4.5) for the blackbody optical 
depth, deviates from numerical solution by less than 1%. 

A very simple correction to the normalization can be arrived at as follows. An immediate 
improvement over the solution of [13] is possible by approximating the non-stationarity in Eq. (3.1) 
using the solution Eq. (3.5) (ignoring the energy injection term) 4 

dn ^ -4fy ^ Cn[(K dC + K br )K c ] l/2 
dt x e 2 dt x e 2 

Note that the time derivative of temperature in Eq. (2.6) can be neglected at small frequencies as its 
effect is suppressed by a factor of x e with respect to the term with the time derivative of yu. Equation 
(3.1) with the above approximation for the non-stationary term gives Bessel's equation 



d 2 d/u (x c 
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The solution is given in terms of modified Bessel function of second kind K v (x), 

fi(x e ) = Afi c J- J—K 05 ^z4c^(x c /x e ). (4.3) 

\ 71 y JCq 

This result provides a more precise dependence of /i on frequency compared to Eq. (3.2). Choosing 
normalization to give /i(0.5) « /u c at x e = 0.5 gives A = 1.007 + 3.5x c , where this fit is accurate 
for 5 x 10~ 3 < x c < 2 x 10~ 2 . This fit thus covers all the interesting range for critical frequency x c 



4 The energy injection term will add a inhomogeneous term to the homogeneous equation for ji(x). The effect of this 
term is to change the overall normalization /j c without significantly affecting the shape of the spectrum. This term can 
therefore be neglected for the purpose of calculating the photon creation/absorption. 
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(see Fig. 3). This choice of normalization frequency (x e = 0.5) provides a good fit to the numerical 
solution. 

The normalization frequency is chosen so that (i) it is < 1 , since this is the assumption made in 
deriving the analytic solution and (ii) it is also large enough so that fj. ~ constant. Since we only want 
to use this solution to calculate the total photon emission, it need only be accurate at x e < 1. The only 
requirement at x e > 1 is that its contribution to the photon emission/absorption should be negligible 
at x e » 1 and that it should be approximately constant around x e = 1, as expected from a correct 
solution. We should also point out that /u(x e ) decreases with increasing x e at x e » 1 for the solution in 
Eq. (4.3) andju(x e — > oo) - 0. Thus our solution satisfies the requirements outlined above. Obviously 
it cannot be normalized at x e = oo, as was done with the original solution Eq. (3.2). The normalization 
must be done by comparison with the numerical solution, taking into account the assumptions made 
in arriving at this solution, resulting in our choice of x e = 0.5. 5 We show a snapshot of the numerical 
solution (chosen at random) at z = 3.48 x 10 6 , original solution Eq. (3.2) and improved solution 
Eq. (4.3) in Fig 4. The critical frequency at this redshift is x c = 0.01 14. Needless to say that the 
shape of the spectrum is well described by our solution at all redshifts and we have chosen a random 
snapshot in Fig. 4 for illustration. The final justification for all our assumptions and approximations 
is of course given by a comparison of the final numerical and analytic solutions for the evolution of 
the chemical potential with redshift as described below. 

We can now use our improved solution to calculate the photon production rate 
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Proceeding as before we get improved formula for blackbody optical depth, 
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3 Since ^(x e ) is approximately constant around x e = 0.5 (variation in the analytic solution is less than 1% for 0.4 < x e < 
1), the exact value of normalization frequency is not important, and we get similar precision if we choose to normalize at a 
slightly different frequency around x c = 0.5, for example at x e = 0.6. We should also mention that the numerical solution 
also shows a tiny decrease (~ 1% from x c = 1 to x e = 100) in the chemical potential at x c > 1 because of the increasing 
efficiency of the recoil effect at high frequencies. 
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Figure 4. Snapshot of the numerical solution, original solution Eq. (3.2) marked SZ1970 and improved 
solution Eq. (4.3) at z - 3.48 x 10 6 (chosen randomly) with the initial // = 1CT 5 at z - 5 x 10 6 including 
both double Compton and bremsstrahlung photon production and no additional energy injection. The critical 
frequency at z - 3.48 x 10 is x c - 0.0114. Both the analytical solutions are plotted with the high frequency 
distortion jj c = 1.39 x 10~ 9 . The numerical and improved solutions are indistinguishable in the figure while the 
original solution underestimates fi{x e ) at low frequencies. 



The improved solution for evolution of fi is still given by the original equation (3.6) on substituting 
the improved optical depth given by Eq. (4.5). 

The improved solution has a broad region of validity and covers the entire redshift range of 
interest. The physics used to derive Eq. (4.5) is applicable for redshifts z < 8 x 10 7 . Thus we can 
use the blackbody optical depth, Eq. (4.5), for redshift interval 10 5 < z < 8 x 10 7 , the upper limit 
is well behind blackbody surface at z ~ 2 x 10 6 . At z ~ 8 x 10 7 the number density of positrons 
becomes comparable to the number density of electrons/baryons due to pair production. At higher 
redshifts, the number density of electrons and positrons, and thus the rates of Compton scattering, 
double Compton and electron-electron and electron-positron bremsstrahlung, increase exponentially 
with increasing redshift. Thus the blackbody optical depth, 7" also starts increasing exponentially 
instead of a power law as in our solution and Eq. (4.5) is no longer applicable. However, we already 
have T ~ 10 4 at z ~ 8 x 10 7 , and creation of a distortion in photon spectrum is thus impossible at 
higher redshifts. 

We plot the optical depth for the numerical solution, the total improved solution as well as 
the individual terms in Eq. (4.5) in Fig. 5. At high redshifts, double Compton terms dominate 
with the new double Compton term Zdc' a l so contributing. Bremsstrahlung term becomes important 
at low redshifts. At high redshifts, our improved solution is indistinguishable from the numerical 
solution. As blackbody optical depth becomes small the quasi-static assumptions made in arriving at 
the analytic solution also breakdown and the error grows. The new solution is however an excellent 
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Figure 5. Blackbody optical depth as a function of redshift (independent of frequency) calculated using numer- 
ical solution, and improved solution Eq. (4.5). Individual terms in Eq. (4.5) are also shown. At high redshifts 
double Compton terms dominate with the new double Compton term Zdc' also contributing. Bremsstrahlung 
term becomes important at low redshifts. At high redshifts our improved solution is indistinguishable from 
the numerical solution. The new solution is an excellent approximation to the numerical result over the entire 
redshift of interest where the optical depth is greater than a few %. 

approximation to the numerical result, and definitely better than the double Compton only result, over 
the entire redshift range of interest, where the optical depth is greater than a few %. The deviations 
from the numerical result for the double Compton only formula and our new result are plotted in Fig. 
6. The analytic solution overestimates the photon production at low redshifts. The reason becomes 
clear by looking at xh in Fig. 3. Photon production at x e > xh would be suppressed since the photon 
production rate is smaller than the expansion rate. At high redshifts, xh > 1 and the error introduced 
by including x e > xh is negligible since photon production is negligible at these frequencies anyway. 
At low redshifts, xh becomes less than unity and starts approaching x c and Eqs. (3.4) and (4.4) 
overestimate photon production. We show the blackbody visibility factors 6 Q = e~ T in Figs. 7 and 8 
for analytic solution Eq. (4.5), for the double Compton only term Zdc and the numerical result. The 
accuracy of blackbody visibility is better than 1% with the new solution Eq. (4.5). 

5 Examples from standard cosmology 

5.1 Upper limit to energy release after BBN and before recombination 

In standard model of cosmology, we can get constraints on energy density in radiation from two 
distinct and very precise observables. The first is the deuterium abundance, which gives the baryon 



6 This is really the visibility of distortions. When the visibility is small, the distortions are not visible, and when the 
visibility is unity, distortions survive and are visible today. 
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Figure 6. Deviation (%) from the numerical solution in blackbody optical depth with respect to the numerical 
solution for the standard double Compton only analytic solution, Eq. (3.7), which includes both bremsstrahlung 
and double Compton using the method of [13] and our new solution. The error in T at low redshifts does not 
have a significant effect on the final spectrum. The error in visibility for our improved solution is better than 
1% at all redshifts and shown in Fig. 8. 

number to photon number ratio r\ = (5.7 + 0.3) x 1(T 10 during primordial nucleosynthesis at 4 x 10 8 > 
z > 4x 10 7 [5, 25]. The second is the measurement of CMB anisotropics, which constrains the baryon 
to photon ratio tj = (6. 1 8 ± 0. 15) x 10~ 10 [6] during recombination at z ~ 1 100. The fact that these two 
values of baryon to photon ratio are almost identical, with small error bars means that we do not have 
arbitrary freedom in adding energy to CMB, for example, with the introduction of new physics. In 
fact any addition of energy/entropy to CMB between primordial nucleosynthesis and recombination 
cannot be more than a small percentage (~ 7% for CMB and BBN to be consistent within 2 - cr) of 
the already existing radiation energy density. COBE limit [24] of fi < 9 x 10~ 5 implies that a ~ 7% 
energy can be added only at redshifts z > 4.1 x 10 6 . These limits also justify our assumption of 
small distortions in the analytic calculations. Any energy injection into photons from non-standard 
processes before electron-positron annihilation, however, is unconstrained. 

5.2 Electron-positron annihilation 

It is, of course, possible to add of order unity energy to radiation before primordial nucleosynthesis. 
This happens in standard cosmology during electron-positron annihilation [32-34], which more than 
doubles the energy density of photons and increases their temperature by ~ 40%. In the early stages 
of electron-positron annihilation, the electrons/positrons far outnumber the photons. In this era, there- 
fore, electron-positron annihilation and electron-electron/electron-positron bremsstrahlung dominate 
the thermalization process. In the very late stages, when most of the electron-positrons have annihi- 
lated, Compton and double Compton scattering are dominant, electron number is conserved and their 
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Figure 7. The characteristics of the blackbody photosphere as given by the blackbody visibility @ = e~ T 
for analytic solution Eq. (4.5), for the double Compton only solution Zdc and the numerical result. We have 
introduced blackbody surface as the boundary where the blackbody optical depth 7~=1. 



density evolves according to the non-relativistic adiabatic law, and our analytic formulae become ap- 
plicable. We, of course, do not expect any observable /u distortion from election positron annihilation 
[13]. It is still interesting to calculate the magnitude of the distortion to demonstrate the effectiveness 
of double Compton scattering and comptonization in restoring the equilibrium between matter and 
radiation. 

At redshifts z < 10 8 , most of the positrons have annihilated and their number density falls 
below that of electrons. The number density of electrons is high enough (as a result of ~ 10~ 9 
asymmetry in matter anti-matter) to maintain the annihilation rate much faster than the expansion 
rate. We can, thus, use Saha equation to follow the positron number density during the last stages 
of positron annihilation. Using the fact that the positron number density is much smaller than the 
electron number density and that the electron number density is unaffected by annihilation, we get 
for the positron number density n+, 

„ 2 



'eq 



where the equilibrium (zero chemical potential) number density of electrons/positrons is 

2 in "'e c2 
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The rate of energy injection is given by, 
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Figure 8. Same as Fig. 7 but going to higher redshifts and also showing the errors, relative to the numerical 
solution, for the different analytic solutions. Deviations from the numerical solution is also shown for Eq. 
(3.7) based on the method of [13] but including both the double Compton and bremsstrahlung. The errors are 
negative at high redshifts and positive at low redshifts with a spike where they change sign. 



We have plotted the resulting /i injection rate multiplied by redshift, BC{\+z) 



2»i e c- dn + „— 7~ 



in Fig. 



9. Visibility function suppresses the high redshift contribution, while the exponentially decreasing 
positron number density suppresses the low redshift contribution, giving the peak at z ~ 1.3 X 10 7 . 
The chemical potential from electron-positron annihilation is suppressed by an astronomical factor 
of 10 178 ! Thus it is impossible to create a deviation from blackbody spectrum at high redshifts. 

5.3 Primordial nucleosynthesis 

Big bang nucleosynthesis (BBN) at z ~ 3 x 10 8 results in binding of almost all neutrons into helium 
(n4 He /nB ~6x 1CT 2 ) along with the production of small amount of deuterium (m H /nB ~ 2 x 10 ), 
helium-3 m He /«B ~ 8xl0 -6 , tritium m H /«B ~ 7xl0 _8 ,beryllium-7/lithium-7 0z7 Be /«B « 3xl0~ 10 ), 
lithium-6 («6 L i/«B * 10~ 14 ) and trace amounts of heavier elements [see 35, for a recent calculation]. 
We can get a rough estimate of the energy released during this main part of nucleosynthesis by 
calculating the total binding energy of helium-4 produced. 7 Thus we have an energy release of 
AE/E ~ Ebi n( im Ue /Ey ~ 6 x 10~ 9 . The blackbody optical depth at z ~ 10 8 is ~ 10 5 and we have the 
final n ~ 0. However, tritium and beryllium-7 survive for a long time before decaying into helium-3 
and lithium-7 respectively. Although, the energy released in the decay of beryllium-7 and lithium-7 
is much smaller than that released during helium production in BBN, the distortions are much larger 
because these decays happen in front of the blackbody surface, when the blackbody visibility is 
almost unity. Also, the energy density released is proportional to the number density of beryllium-7 



7 We ignore the fact that some of the energy will be lost to neutrinos. 
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Figure 9. Chemical potential fi from electron positron annihilation. The CMB blackbody spectrum is main- 
tained at extraordinary precision of 10~ 178 ! Also shown, in the inset, is the actual rate of energy injection 
multiplied by (1 + z). At high redshifts (z > 10 8 ) we have used entropy conservation to calculate the rate of 
heating. 



and lithium-7, which has decreased as (1 + z) 3 compared to the (1 + z) 4 decrease of the radiation 
energy density, thus giving a larger AE/E than if the decay had happened at the same time as the 
main BBN. Tritium has a half life of 12.32 years. It, therefore, decays at z ~ 2.5 x 10 5 to helium-3 
releasing an electron with average energy 5.7KeV. Most of this energy release happens at z > 10 5 and 
causes a //-distortion with/i = 2 x 10 -15 . 

Neutral beryllium atom decays by electron capture with a half life of 53.2 days. Fully ionized 
beryllium in the low density plasma in the early Universe is however stable. It has to wait for « 800 
years until z ~ 3 x 10 when it recombines to hydrogen like beryllium. The recombined beryllium 
can now capture the orbital electron and decay to lithium-7 with a half-life of 106.4 days, which 
is twice the half-life of a fully recombined beryllium [36]. 89.6% of the decays go to the ground 
state of lithium-7 and most of the energy released is carried away by neutrinos, which would appear 
today as a narrow line in the cosmic neutrino spectrum. 10.4% of beryllium decays into an excited 
state of lithium. About half of the total decay energy in this case also is lost to neutrinos forming 
a second lower energy line in the cosmic neutrino spectrum. The excited lithium nucleus then de- 
excites, almost immediately, to the ground state, emitting a Q = 477.6 KeV photon, which delivers 
most of its energy to plasma by Compton scattering on electrons (recoil effect). The Compton y 
parameter at z = 30000 is 0.04, which lies intermediate between pure y and jj. type eras. The heating, 
therefore, results in a distortion intermediate between the y and the /j type distortions of magnitude 
(using formula for j-type distortion) « (l/4)AE/E ~ (l/4)0.1042«7 Be /£' y ~ 10~ 16 . 



5.4 Dark matter annihilation 

A natural and favored candidate for dark matter is a weakly interacting massive particle (WIMP) with 
several candidates in high energy theories beyond the standard model [37]. A very attractive feature 
of WIMP is that if they have weak scale interactions then the correct amount of dark matter (which is 
close to the critical density) observed today can be thermally produced in the early Universe. This is 
remarkable since a priori there is no reason to suspect any relation between the weak scale interactions 
and the present critical density of the Universe and this coincidence is sometimes referred to as the 
WIMP miracle. For the thermally produced WIMPs, the dark matter density is related to the velocity 
averaged cross section as 

, , 3x 1(T 27 3 .J 

<crv) « — em's , (5.4) 

O dm /^ 

where ho « 0.702 is the Hubble parameter and Q am is the dark matter density as a fraction of critical 
density today. The above values for annihilation cross sections are of similar order of magnitude 
as the current upper limits from Fermi-LAT experiment probing dark matter annihilation in the lo- 
cal Universe [38^4-0]. The actual annihilation cross section can, of course, be much smaller than 
the Fermi upper limits. Energy released from dark matter annihilation also changes recombination 
history [41] and it can thus be constrained through its effect on the CMB power spectrum [42, 43] 
and recombination spectrum [44]. Effect of dark matter annihilation also changes the abundance 
of elements produced during BBN [45, 46]; these constraints are complementary but less stringent 
compared to the current CMB constraints. 

Initially, dark matter is in thermal equilibrium with other constituents of the Universe, and is 
being continuously created and annihilated. As the Universe cools and expands, these interactions 
freeze out and thereafter the dark matter number is conserved. However, a small number of residual 
annihilations keep happening throughout the histoiy of the Universe. The rate of energy released into 
CMB from these residual annihilations 8 is given by, 
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where tmwimp is th e mass of the dark matter particle and f y is the fraction of energy that goes into 
particles with electromagnetic interactions, and is deposited in the plasma. 
We plot the rate of energy release into the CMB 

(l+z)g^=S e — (5.6) 

dz H 

in Fig. 10 for the fraction of energy going into the plasma 9 f y = 1. Hubble rate is proportional to 
(1 + z) 2 during radiation domination, which gives us the flat portion of the curve. Energy release 
rate decreases faster than the expansion rate <x (1 + z) 3/2 during matter domination, giving the low 
redshift declining tail in the plot. The total ji distortion (using energy release from z > 5 x 10 4 ) is 



We assume self annihilating Majorana particles. For Dirac particles the energy release would be smaller by a factor of 
2. 

'in general we expect some energy to be lost to neutrinos and other dark particles and therefore f y would be less than 
unity. 
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Figure 10. Energy injection from dark matter annihilation for a 10 GeV WIMP with f y — 1 (solid line). The 
exponential suppression at high redshifts is because of the decrease in visibility, e~ T . Dashed line shows the 
energy injection (1 + z)% without the visibility factor. 

fi v 3 x 1(T 9 . The y- type distortion from energy release z < 5 x 10 4 is y « 5 x 10 -10 , these distortions 
were also calculated by [23]. These numbers are of similar order of magnitude as the distortions from 
Silk damping and Bose-Einstein condensation of CMB discussed in the next two sections. COBE 
constraint of /i < 9 x 1(T 5 [24] constrains WIMP mass to be m W i M p > 0.3/ y MeV while PIXIE [27] 
would be able to constrain up to wiwimp ~ 3/ r GeV. 

5.5 Silk damping 

Sound waves are excited in the primordial baryon-electron-photon plasma by primordial perturba- 
tions. They decay on small scales because of shear viscosity, with thermal conduction also becoming 
important near the time of recombination. This damping of primordial perturbations was first cal- 
culated by Silk [47] including only thermal conduction. The full calculation, including both shear 
viscosity and thermal conduction and also including the effects of photon polarization was done by 
Kaiser [48]. Silk damping transfers energy from sound waves to the average CMB spectrum, result- 
ing in effective energy injection into CMB [23, 49-52]. Microscopically, shear viscosity and thermal 
conduction arise due to the diffusion of photons which are repeatedly scattered by the electrons. This 
diffusion of photon results in mixing of blackbodies from different phases of the sound waves on 
diffusion scales. 2/3 of the dissipated energy in sound waves just increases the average temperature 
of CMB while 1/3 results in the spectral distortions of yu and y type. Depending on the primordial 
perturbation power spectrum at these very small scales of comoving wavenumbers 50 < k < 10 4 , the 
fi distortion can be in the range 10~ 8 - 10~ 10 [52] for the parameter space allowed in the standard 
cosmological model [6]. We refer to [52] for a detailed discussion, including fitting formulae for 
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Process 


H 


electron-positron annihilation 


10 -178 


BBN tritium decay 


2 x 1(T 15 


BBN 7 Be decay 


1Q -16 


WIMP dark matter annihilation 


3xlO- 9 /y 10GeV 


Silk damping 


10-8 _ 10 -iJ 


Adiabatic cooling of matter and 




Bose-Einstein condensation 


-2.7 x 10" 9 



Table 1. Census of energy release and fi distortions in standard cosmological model. The negative distortion 
from adiabatic cooling of matter is shown in red. 



Process 


y 


WIMP dark matter annihilation 


6xl(r 10 /yL 0UeV 


Silk damping 


ifr 8 - io-» 


Adiabatic cooling of matter and 




Bose-Einstein condensation 


-6 x ifr 10 


Reionization 


lO" 7 


Mixing of blackbodies: CMB € > 2 multipoles 


8 x ifr 10 



Table 2. Census of energy release and y distortions in standard cosmological model. We also give the value 
of y-type distortion expected from the mixing of blackbodies when averaging our CMB sky [53]. The negative 
distortion from adiabatic cooling of matter is shown in red. y type distortion is clearly dominated by the 
contributions, during and after reionization, from the intergalactic medium and clusters of galaxies, and the 
early Universe contributions are difficult to constrain. 

spectral distortions from adiabatic initial conditions, and constraints from the future experiments on 
initial power spectrum spectral index and its running. 

5.6 Bose-Einstein condensation of CMB 

After the epoch of electron-positron annihilation, electrons and baryons are non-relativistic and cool 
adiabatically (with adiabatic index 5/3) as a result of the expansion of the Universe, T e oc (1 + z) 2 . 
Radiation (photons) has adiabatic index 4/3 and cools slower than baryons, T y oc (1 + z) [2], Comp- 
tonization however is very efficient before recombination and efficiently transfer energy from pho- 
tons to electrons/baryons, keeping them at same temperature as photons. This cooling of CMB 
[23], along with thermalization from comptonization, results in Bose-Einstein condensation of CMB 
[26]. The photons thus move from high to low frequencies where they are efficiently destroyed by 
bremsstrahlung (and at high redshifts also by double Compton scattering). Since the amount of cool- 
ing is small, linear theory for small distortions applies. The resulting distortions have the same shape 
as that caused by heating of CMB in previous examples, but with opposite sign. Thus we have nega- 
tive n and negative y distortions which partially cancel the distortions due to dark matter annihilation 
and Silk damping. Surprisingly, the // (and y) distortions have a magnitude which is similar to those 
from dark matter annihilation and Silk damping. A comparison of /i distortions from Bose-Einstein 
condensation as well as all previous examples is presented in Table 1. We also show comparison 
of y-type distortions in Table 2. y-type distortions are dominated by the low redshift contributions, 
during and after reionization, from the intergalactic medium and clusters. Early universe physics is 
therefore difficult to constrain using the y-type distortions. 
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5.7 Energy released from recombination of plasma 

We should also mention that recombination lines also create a distortion of amplitude AT IT ~ 1(T 8 - 
10 -9 [see 54, for a detailed calculation]. The distorted spectrum would heat the electrons, adding a 
y-distortion at the time of Helll— >HeII recombination of ~ y(6000) x 10~ 9 ~ 10~ 12 . Additionally, 
the Ly-a and 2s- Is (2-photon decay) photons from recombination with energy ~ 40eV, x ~ 30, 
escape as they redshift out of resonance, (Compton) scatter on electrons and heat the plasma through 
recoil effect. The heating can be estimated using the analytic solution of Kompaneets equation with 
only the recoil effect [55, 56]. In the limit of small Compton-y parameter (xy «: 1), the fraction of 
energy lost by photons at frequency x is ~ y x x ~ 1/30, giving an additional y-type distortion of 
~ (l/4)(40eV)m He xy/£' y ~ 10~ 12 . The distortions from Hel and HI recombination are much smaller 
since they happen later, when y is much smaller, although the energy released is comparable to Hell 
recombination. 



6 Conclusions 

Future experiments, such as PIXIE [27], would be able to constrain/measure spectral distortions in 
the CMB at high accuracy. There are several sources of spectral distortions possible from standard 
and new physics. Using the results of experiments like PIXIE to constrain new physics would 
require precision calculations of evolution of the CMB spectrum, especially around the blackbody 
surface at z ~ 2 x 10 6 . So far, precise calculations have only been possible numerically, although 
analytic solutions with 5 - 10% precision around the blackbody surface have been available for a 
long time. We have presented new analytic solutions, which take into account both double Compton 
scattering (important at high redshifts) and bremsstrahlung (important at low redshifts). We also 
take into account the non-stationarity of the problem which is important to achieve high precision. 
The new solutions are presented in Eq. (3.7) (ignoring non-stationarity) and in Eq. (4.5) (including 
the non-stationarity of the problem). Equation (4.5) gives accuracy of better than 1% in blackbody 
visibility at all redshifts. We also present examples from the standard AC DM model of cosmology, 
which do not require new physics, illustrating the structure of blackbody surface. In particular, 
electron-positron annihilation and BBN demonstrate the effectiveness of Compton and double 
Compton scattering in maintaining equilibrium at high redshifts. We also point out the 
coincidence/degeneracy among the three significant sources of distortions in standard cosmology, 
Bose-Einstein condensation of CMB, Silk damping and dark matter annihilation. All of these create 
distortions which have roughly the same order of magnitude, especially for low dark matter particle 
masses. This is remarkable considering that the three sources of distortions have completely 
different physical origins. Bose-Einstein condensation is fixed by standard cosmological parameters, 
which are now known with high precision. However, the degeneracy between Silk damping and dark 
matter annihilation must be taken into account when using spectral distortions to constrain the 
primordial power spectrum or dark matter parameters. 
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